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?3 ■ Abstract 

We consider the weak solution of the simphfied Ericksen-Leslie system mod- 
ehng compressible nematic liquid crystal flows in M^. When the initial data is 
^> \ small in and initial density is positive and essentially bounded, we first prove 

■ the existence of a global weak solution in M'^. The large-time behavior of a global 

^ . weak solution is also established. 

d 

(N 

1 Introduction 

^ ! We consider the following hydrodynamic system modeling the flow of nematic liquid 

ci\ crystal materials [2,7,20], 

Pt + V-(pn) = 0, (1.1a) 

put + pu-Vu + VP{p) = pAu + XVdivu - Vd ■ Ad, (1.16) 

dtd + u-Vd = Ad+\Vd\^d, (1.1c) 

for {t,x) G [0,+oo) X M^. Here p, u = {u^ ,v? ,v?y and P denote the density, the 
velocity, and the pressure respectively, d = {d^,d'^,d^y is the unit-vector (\d\ = 1) on 
the sphere §^ C representing the macroscopic molecular orientation of the liquid 
crystal materials, p and A are positive viscosity constants, and div and A are the usual 
spatial divergence and Laplace operators. 

The above system (1.1) is a simplified version of the Ericksen-Leslie model for the 
hydrodynamics of nematic liquid crystals. The mathematical analysis of the incom- 
pressible liquid crystal flows was initialed by Lin and Liu in [21,22]. In [29,30], Wang 
proved the global existence of strong solutions in whole space under some small con- 
ditions. In dimension two, the global existence and uniqueness of weak solutions were 
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studied in [10,11,18,23,31] and reference therein. On the other hand, the global exis- 
tence of weak solutions to incompressible liquid crystal in is still an outstanding 
open question. 

When the fluid is allowed to be compressible, the Ericksen- Leslie system becomes 
more complicate. To our knowledge, there seems very few analytic works available 
yet. The local-in-time strong solutions to the initial value or initial boundary value 
problem of system (1.1) with nonnegative initial density were studied in [3,13]. Based 
on [16,17], the blow up criterion of strong solutions were obtained in [13,14]. The 
global existence and uniqueness of strong solution in critical space were studied in [12]. 
Motivated by [15], when the initial data were sufficiently smooth and were suitably 
small in some energy-norm, the global well-posedness of classical solutions were proved 
in [19]. Especially global weak solutions in two dimensions was established in [6] under 
some small condition or geometric angle condition. 

Our aim in this paper is to establish the global existence of low-energy weak solu- 
tions of system (1.1), with the following initial conditions: 

(p(-,0),«(-,0),d(-,0)) = (po,wo,c?o), (1-2) 

where po is bounded above and below away from zero, \do\ — 1, Mq, V(io £ L^(]R^) for 
some p > 6, and modulo constants, (po, Uo, do) is small in L^(M^). Thus the total initial 
energy is small, but no other smallness or regularity conditions are imposed. 

When the direction field d does not appear, (1.1) reduces to the compressible Navier- 
Stokes equations. The global classical solutions were first obtained by Matsumura- 
Nishida [25,26] for initial data close to a non- vacuum equihbrium in H^(R^). In par- 
ticular, the theory rcqiiircs that the solution has small oscillations from a uniform 
non-vacuum state so that the density is strictly away from the vacuum and the gradi- 
ent of the density remains bounded uniformly in time. Later , Hoff [8,9] studied the 
problem for discontinuous initial data. For the existence of solutions for arbitrary data, 
the major breakthrough is due to Lions [24] (see also Feireisl [5]), where he obtains 
global existence of weak solutions-defined as solutions with finite energy. Suen and Hoff 
[28] adopted Hoff's techniques to obtain global existence of low-energy weak solutions 
for the magnetohydrodynamics. In this paper, we shall study the Cauchy problem 
(1.1)-(1.2) for hquid crystals and estabhsh the global existence and large time behavior 
of low-energy weak solutions However, compared with the compressible Navier-Stokes 
equations, some new difficulties arise due to the additional presence of the liquid crystal 
directional field. Especially, the super critical nonlincarity |V(ip(i in the transported 
heat flow of harmonic map equation (1.1c) and the strong coupling nonlinearity Vd-Ad 
in the momentum equations (1.16) will cause serious difficulties in the proofs of the 
time-independent global energy estimates. 

To state the main results in a precise way, we first introduce some notations and 
conventions which will be used throughout the paper. For a given unit vector n e 
and a positive integer m, we denote 

^m(^3. g^).-{d:d-ne H'^iR^), \d\ = 1 a.e. in M^}. 
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We use the usual notation for Holder seminorms: for v : — > M.^ and a e (0, 1], 

\V{X2) - V[xi)\ 



{v)"' — sup 



,X2 tK"- 



\X2 — X\ 



a 





\V{X2, t 


2) - v{xi,ti) 


1 


1X2 - Xi 


\ai _|_ 


^2 — ^1 


0-2 



and for : g C M3 X [Q, 00) — > and ai, cts e (0, 1], 
{v)q = sup 

If X is a Banach space we will abbreviate by X when convenient. Finally if 
/ C [0,00) is an interval, C^(/;X) will be the elements v e C(/;X) such that the 
distribution derivative Vt G V'{int /; M^) is realized as an element of C(/; X). 

As it was pointed out in [8], the effective viscous flux plays an important role in 
the mathematical theory of compressible fluid dynamics. More precisely, let F and uj 
be the effective flux and vorticity deflned by 

F = (/X + \)divu - {P{p) - P{p)) and cu = V xu. (1.3) 

It is not hard to check that 

So, it follows from (1.1b) that 

AF = (iz't;(pM + rfw(Vrf(8) Vrf) - V^^^), //Au; = V x (pu + (iii;(Vc?(8) Vo?)), (1.4) 
where " ' " denotes the material derivative, i.e.. 

Now we give a precise formulation of our results. First concerning the pressure P, 
we focus our interest on the case of isentropic flows and assume that 

P(p)=ap^ with a>0,7>l. (1.5) 

Next we flx a positive reference density p and then choose positive bounding densities 
p and p satisfying 

p < min{p, p'} and max{p, p"} < p, (1.6) 
and flnally we deflne a positive number 5 by 

S = min{min{p, p'} - p,p- max{p, p"}, ^{p - p)}. (1.7) 



(Notice that S need not be "small" in the usual sense.) Concerning the diffusion 
coefficients /i and A we assume that 

< A < ^±^;,. (1.8) 

It follows that 

for p = 6 and consequently for some p > Q, which we now fix. 

Concerning the initial data (po,^0)'^o)i we assume there is a positive number TV, 
which may be arbitrary large, and a positive number b < 6 such that 

\\uo\\Lp + \\Vdo\\LP < N (1.10) 

and 

p + b < essinf po < ess sup po < p — b. (l-H) 

We assume also that 

do-ne H^{R^;8^) (1.12) 

and write 

Co= / {lpo\uo\' + G{po) + l-\Vdo\'')dx, (1.13) 
where G{p) is the potential energy density defined by 

G(p) = p£Ml^.. (1,4) 

It is clear that there exist two positive constant ci,C2 only depending on p,p, and p 

c,{p,p,p){p-pY < G{p) < C2ip,p,p){p-p)'. (1.15) 

Weak solutions of (1.1)-(1.2) are defined in a usual way. 

Definition 1.1. A pair of functions {p,u,H) is said to be a weak solution of (1.1)- 
(1.2) provided that {p - p, pu) E C{[0, oo); H-\R^)), d-n e C{[0,oo); L'^{m.^)), 
u,Vd e L°°([0,oo);L2(M3))^ Vu G L'^{{0, oo); L^{R^)), and \d{-,t)\ = 1 a.e. in for 
t >0. Moreover, the following identities hold for any test function ip e P(R^ x (^1,^2)) 
with ^2 > ^1 > and j = 1, 2, 3: 

/ p^{x,t)dx\ll — / / {pipt + pu ■ 'Vtp)dxdt, 

Jr3 Jt^ J^3 

/r3 pu^ij(x,t)dx\ll + /jg3(/^Vii^ • VV' + X(divu)'ijja:j)dxdt 
= It' /m3 (P^'i't + pu^u ■ Vt/j + P{p)^^, - \ I Vd\^j + 4, Vd • VV') dxdt, 

4 



and 

J^aid^ — n^){x, t)'i/j{x, t)dx\ll + It ^ /r3 ^d^^ipdxdt 
= JT*' J^s{{d^ - n^)i)t - u ■ VdH) + \vd\'^d^il:)dxdt. 
Our main results are formulated as the following theorem. 

Theorem 1.2. Assume that the system parameters in (1.1) satisfy the conditions 
(1.5)-(1.9) and let positive numbers N and h < 5 he given. Then there are positive 
constants e, C, and 9 depending on the parameters and assumptions in (1.5)-(1.9), on 
A^, and on a positive lower bound for 6, such that, if initial data [pQ^UQ^do) are given 
satisfying (1.10)-(1.14) with 

Co < £, (1.16) 

then there is a solution {p,u,d) to (1.1)-(1.2) in the sense of the definition 1.1. More- 
over, the solution satisfies the following: 

p-p,pueC{[Q,ooy,H-\R')), (1.17) 

d-neC([0,oo);L\M.^)), (1.18) 

Vm e L2([0,oo);R^), (1.19) 

u(-,t),Vd(;t) e H\R^), t>0 (1.20) 

F(;t),co(;t) e H\R^), t>0 (1.21) 

(«)M[,oo)> (Vc^)M[,^) < C(r)Co^ (1.22) 
where C(t) may depend additionally on a positive lower bound for r, 

p < p{x, t) <p a.e. on x [0, oo), (1.23) 

and 

sup J^,[\p - + \u\^ + i|Vdp + ai\Vu\^ + iV^dp) + a5(F2 + \Vuf)]dx 



TOO 



+ r /m3[|V«P + I Ad + \Vdm' + a{\u\' + |Vd,p + iVcp) n .24) 

+a5(|VM|2+ |V2(ii|2)](ix(is 
where a{t) = min{l,t}. Moreover we also have the following large-time behavior: 

hm (Hp - p||i,((]R3) + ||w||vyl,r(]K3) + ||Vd||^l,r(K3)) = 0, (1.25) 

holds for ; e (2,oo),r e (2,6). 

The rest of this paper is devoted to prove Theorem 1.2. In Section 2, we collect 
some useful inequalities and basic results. In Section 3, we derive the time-independent 
energy estimates of the solution. The key pointwise upper and lower bound of the 
density are established in Section 4. Finally, the proof of the main results will be done 
in Section 5. 
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2 Preliminaries 

In this section, we state some auxiliary lemmas, which will be frequently used in the 
sequel. We start with the well-known Gagliardo-Nirenberg inequality (see, for instance, 
[1,32]). 

Lemma 2.1. First, given r e [2, 6] there is a constant C(r) such that for / e 

II/I|l^(m3) < C{r)\\f\\%-^J^'^\\Vf\\^^^^^^^^^ (2.1) 
Next, for any r G (3, oo) and q > 1, there is a constant C{r,q) such that for / G 

||/IUoc(«3) < CM\\f\\t[^J^/^^^^^^^^ (2.2) 

and 

(/)R3<C(r)||V/|U.(M3), (2.3) 

where a = 1 — -. 

r 

The next lemma is due to Hoff [8], which will be used to prove the uniform (in 
time) bound of density. 

Lemma 2.2. If F is the fundamental solution for the Laplace operator in M^, then 
given pi e [1, 3) and p2 G (3, oo], there is a constant C — C{n,pi,p2) such that 

llr^,- * /||l-(M3) < C(n,pi,p2)[\\f\\LPi{R3) + ||/||lp2(R3)]. (2.4) 



Finally, we need the local-in-time existence theorem of (1.1)-(1.2). Using the mol- 
lifier technique, the local solutions which can be proved rigorously by the standard 
method of Matsumura and Nishida [25,26]. 

Proposition 2.3. Assume that the initial data {po,Uo,do) satisfies 

(po - p, uo) e H%W^) do e HW; §2) and inf (po(x)) > 0. (2.5) 

a;6R3 

Then there exists a positive time To, which may depend on mf^^^si^po^x)), such that 

the Cauchy problem (1.1)-(1.2) has a unique smooth solution {p,u,d) on R? x [0,71)] 
satisfying 

p{x,t) > for all X G M^t e [0,To], (2.6) 

p - p e C([0, To]; i/3) n C\[Q, To];H'), (2.7) 

u e C([0, To]; H') n C\[0, To]; H'), (2.8) 
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d-ne C([0, To]; H^{R'; S^)) n C\[0, To]; i7^(R^ S^)). (2.9) 
In view of Lemma 2.1 and the classical estimates of elliptic system, we have 

Lemma 2.4. Let {p, u, d) be as in Proposition 2.3, then there exists a generic positive 
constant C, depending only on /i and A and p, such that for p e (1, oo) 

||VT||ip(]R3) + ||Va;||iP(K3) < C(||pm||^p(]r3) + ||Vc? • Ad\\LP(j^3)), (2.10) 

||Vm||z,p(k3) < C(||T||^p(K3) + ||a;||LP(K3) + \\P{p) - T'(p)||lp(r3)), (2.11) 
where F and u are defined in (1.4). 

Proof. An apphcation of the i/-estimate of elhptic systems to (1.4) gives (2.10). 
On the other hand, since —Au — —Vdivu + V x a;, it holds that 

Vm = -V{-A)-^Vdwu + V(-A)- V X w, 

which, combined with the Marcinkiewicz multiplier theorem (Stein [20,p.96]), we arrive 
at 

||Vri||LP(M3) < C{\\divu\\LP{w:i)) + ||^||lp(r3) 

< C(||F||ip(K3) + ||w||lp(r3) + \\P{p) - P(p)||iP(K3)). 

Thus the proof the lemma is completed. 



3 A priori estimates 

This section is devoted to establish a number of a priori bounds for local-in-time 
smooth solutions, corresponding roughly to (1.24). Those are rather long and technical. 
We have therefore omitted those which are identical to or nearly identical to arguments 
given elsewhere in the literature of whose details we regard as routine. On the other 
hand, we have endeavored to describe the flow of the arguments in such a way that the 
diligent reader can reconstruct the details without undue difficulty. 

Let T > be fixed and assume that (p, u, d) is a smooth solution of (1.1)-(1.2). We 
define a functional A{t) for a given such solution that 

A{t) = sup J^s[a{\Vu\^ + iV^dp) + c75(|m|2 + \Vu\^ + \Vdt\^)]dx 

o<s<t (3.1) 

+ /o /M3k(l«l' + I'^dt? + |Vu;n + a%\Vu\^ + \V^dt\^)]dxds, 

where a{t) = min{l, t}, and we obtain the following a priori bound for A[t) under the 
assumptions that the initial energy Cq in (1-13) is small enough and that the density 
remains bounded above and below away from zero: 

Proposition 3.1. Assume that the system parameters in (1.1) satisfy the condi- 
tions in (1.5)-(1.8) and let positive numbers N and b < 5 he given. Assume {p,u,d) 
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is a solution of (1.1) on M.^ x [0,T] in the sense of Proposition 2.3 with initial data 
{po,Uo) e H'^{R^) and (Iq E H^{M.^;S^) satisfying (1.10)-(1.13), then there are positive 
constants e, M, and 6 depending on the parameters and assumptions in (1.5)-(1.8), on 
A^, and a positive lower bound for b, such that if Co < e and 

P < p{x,t) <ponR^x [0,T], 

then 

A{T) < MC^. 

The proof will be given in a sequence of lemmas in which we estimate a number of 
auxiliary functionals. To describe these we first recall the definition of (1.9) of jo, which 
is an open condition, and which therefore allows us to choose q e [6, min{p, 12}) which 
also satisfies (1.9). Then for given {p,u,d) we define 

A(t)^ sup J^s{\^u\'^ +\V^d\'^ + \u\'^ + \Vujf + \Vdt\'^)dx 

+ /o /m3(I^I' + iVrftP + |Va;P + \Vu\^ + \V^dt\'')dxds, 

Bq{t) = sup L{\u\i + \Vd\'})dx + /* L{\u\i-^\Vu\^ + \Vd\i-^\V^d\^)dxds 

0<s<t 

+ io7M3(l«l^-'|V(|«P)P+|Vci|'^-^|V(|Vdp)P)dxd., 
D{t)^ sup f (T5(|VdnV'd|2 + |^/|2|V'd|2 + |VunVd|')dx, 

0<s<t Jm.3 

D{t)^ sup [ (|Vc/HV'd|' + |^/nV'c/|2 + |ViinV(/|')dx, 

l<S<t Jr3 

E{t) = J^sWH\Vu\' + |V2d|3) + a^dV^I^ + \W^^)]dxds 
+1 E IolR^<^<,<,</^ds\, 

and ^ 

E{t)= I I {\Vu\^ + \Vd\^\V^d\^ + \Vu\'' + \V^d\'')dxds. 

It will be seen that the assumed regularity (2.6)-(2.9) suffices to justify the esti- 
mates that follow. We begin with the following energy estimate. 

Lemma 3.2. Assume that the hypotheses and notations of Proposition 3.1 are in 
force. Then 

sup J^adp - p\'^ + \u\^ + \Vd\^)dx 

o<t<T (3.2) 
+ J^J^,{\Vu\^ + \Ad+\Vd\^d\^)dxdt < MCo. 

Proof. Multiplying (1.1b) by u and integrating over R^, we have 
^ / l-p\u^\dx+ f VP{p)udx+ f {fx\Vu\'^ + X\divu\'^)dx ^ - f u-Vd-Ad. (3.3) 
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By the mass equation (1.1a) and the definition of G{p) in (1.14), we have 

G{p)t + div{G{p)u) + (P(p) - P)div u = 0. 
Integrating and adding the result to (3.3) we obtain 

4 / IpIu"^] + G{p)dx + [ {p\Vuf + X\divu\'^)dx^ - f u-Vd-Ad. (3.4) 

dt Jr3 2 J]r3 J-g3 

Multiplying (1.1c) by Ad+ \'Vd\'^d and integrating over M^, using integration by parts 
and the fact that \d\ = 1 we obtain 

T" / hvd\'^dx+ f \Ad+\Vd\'^d\'^dx = f u-Vd-Ad. (3.5) 



dt Jk3 2 

Adding (3.4) to (3.5) and integrating over [0,t], yields (3.2) by (1.15). Thus the proof 
of lemma is completed. 

Lemma 3.3. Assume that the hypotheses and notations of Proposition 3.1 are in 
force. Then for < t < 1 A T 

sup (7 /i83(|Vii|2 ^ \^'id\^)dx + a{\u\^ + \Vdt\^)dxds 

0<s<t (o a.) 

< M[Co + Of' Bf' + C^~' Bl' + E] , 
and if T > 1 and 1 < t < T, then 

sup /jg.,(|Vtip + \V'^d\'^)dx + /o /r3(I^P + \^dt\^)dxds 
i<^<* 3 _ _ (3.7) 
< M[Co + C|A5) + ^] + A(l). 



Proof. For < t < 1 A T, multiplying the equation (1.1b) by aii and integrating over 
X [0,t], we have 

sup a L3 \Vu\'^dx + Jn Ls \u\^dxds 

o<s<t (3.8) 
< M{Co + I /o a[u{div{Vd ® Vrf) - ^V\Vd\^)dxds\ + £;}. 

Differentiating (1.1c) with respect to x, we have 

Vdt- AVd^V{\Vd\'^d-u-Vd). (3.9) 

Multiplying the above equation by aVdt and integrating over x [0,t], we have 

|(j Jjj3 \Ad\^dx + Jjg3 <t| Vcftpc/xds 
= I /o /ig3 o-'l A(i|2(ia;c?s + /J /]g3 (7VdtV(| Vdpd - u ■ Vd)dxds. 



Adding this to (3.8) and combining with Cauchy's inequahty we then get 

sup a /]k3(|Vm|2 ^ I A(i|2)c/a; + /J /^g a{\u\^ + \V dt\'^)dxds 

0<s<t 

< M{Co + E + J^J^,a'\Ad\^dxds (^-lO) 
+ /o Lkd V2(f|2| Vdp + I Vupl Vdp + I V2ci|>|2)](ixd4. 

The right terms can be estimated as follows: 

/o J^s (T'\Ad\^dxds < M Jl""^ /j^3 \Ad + |V(ipd|2 + |V(i|^(ixc/s 

< M[Co + /m3 I VdPdxd.)!^ /J^^ \Vd\'^dxds)^^] 
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<M(Co + Co^-^5r^), 

SlS^s<7\Vu\''\Vd\''dxds < M{jlj^,a^Vu\Hxds)hi'j^, \Vd\'dxds)h 

< M[i/+ (/y^3 \Vd\HxdsY^^{fj^, \Vd\<^dxds)^^] 

< M(Co + c|^S|^). 

The other two terms in the integral on the right side of (3.10) are bounded in a similar 
way, and (3.6) follows. 

For 1 < i < T, as in (3.10), we have 

sup J^,{\Vu\'^ + \Ad\'^)dx + J^,{\u\'^ +\Vdt\^)dxds 

l<s<t 

< M{Co + E + JI /k3[(| V^dH Vc^p + I Vc/|2 + \W^d\^\u\^)]dxds} + A{1). 

(3.11) 

Using the fact that 

I Vrfp = -d-Ad {since \d\ = 1), (3.12) 
the right terms can be bounded as follows: 



Ii Irs \Vu\^\Vd\^dxds < M Jl /jj3 1 V^xp + \Vd\^dxds 

< M [I /„3 1 V^xP + \Vd\''\V''d\''dxds 

< ME 



IiIrs \V^d\^\u\^dxds < \V^d\^ + \u\''dxds 

< ^ + /i(/r3 H^dx)HlM3 \^u\^dx)-2ds 

< m{e + cIa^). 

Taking the above results into (3.11), then (3.7) follows. Thus the proof of lemma is 
completed. 

Next we derive prehminary bounds for ii and Vdt in L°°{[0, T]; L^(M^)). 
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Lemma 3.4. Assume that the hypotheses and notations of Proposition 3.1 are in 
force. Then for < t < 1 A T, 



sup /^3(|u|2 + \Vdt\^)dx + /o J^, (j^(|Vm|2 + \V^dt\^)dxds 

0<s<t 

< M[Cq + E + C'^sf^ + cl=^Bf^ + ct~^Bf~\E + c'f^Bf^) (3-13) 

+(c|^sr)U], 

and if T > 1 and I < t <T, then 

sup jRsditp + \Vdt\^)dx + /o /RsdVitp + \W^dt\'^)dxds 



l<s<t 

2(g-3) 



(3.14) 



< M{Co + CoAE + Co'*^-'' Bg'^'-'^ A + E} + A{1). 



Proof. By the definition of material derivative, we can write (1.2) as follows, 

pii + V(P(p)) = f^Au + Wdivu - Vd ■ Ad. (3.15) 

Differentiation (3.15) with respect to t and using (1.1), we have 

put + pu-Vu + V{P{p)t) + {Vd ■ Ad)t 
= pAu + \Vdivu-[pA{u-Vu) + \Vdiv{u-Vu)] (3.16) 
+div[{pAu + XsJdivu) ®u- VP{p) ®u - {Vd ■ Ad) ® u]. 

Multiplying (3.16) by a^ii, and integrating over x [0,t], we obtain that for < t < 
1 AT, 

sup (7^ Jjja l^pfia; + j^3G^\Vu\'^dxds 

0<s<t 

< M[Co + E + Bg'-' + B^^ ^'^■^'> 
+ /o /r3 ^'|Vrf|2(|Vrft|2 + \u\''\V''d\'')dxds\. 

Next we differentiate (3.9) with respect to t, multiply by a^Vdt and integrate over 
W X [0,t] to obtain 

/r3 \ Vdt\^dx + jl /jg3 a^\V'^dt\^dxds 
= I /o /jj3 a^a'\Vdt\Hxds + /^s a^V(| Vd^c^ - u ■ Vd)tVdtdxds. 

Adding this to (3.17), integrating by parts, using Cauchy's inequality, we then have 



sup a' J^s \u? + \^d\^dx + £ (J^d Vu|2 + \W^d\^)dxds 

0<s<t 

< M[Co + E + Ci^' Bf'' + cl^' Bf^ (3.18) 
+ /o J^, a5|V(^|>p(|VM|2 + iV^rfp) + a^\Vd\yt\^dxds 
+ /o Irs (j'(|Vc/|>|2 + iVdPlVc/tP + \Vdt\'\u\'')dxds]- 
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By (2.2), the terms on right side can be bounded by 



/o /r3 cr^|V(ip|-up|Vnp(ix(is 

< J^J^,a^{\Vd\^ + \u\^ + \Vu\'')dxds 

< E+ sup ||KV(i)G9)||^,(j,3)/oV5||(w,V(i)(s)||i^(j,3)(i5 



0<s<t 



By (1.1c) and (3.12), we have 

< M /J ^3 a^d Vrfh V^dp + I VdHtipl Vd|2 + \Vd\^)dxds 

< M /J jr^3 a^d Vrfh V2rf|2 + \V^d\M^\Vdf)dxds 

< ME + M[C^ Bf^] [E + ] . 
The last term on the right side in (3.18) can be bounded by 

^ (/o t |W|'^cixds)i(/g 43(T¥|M|3cia:cis)i 

< (Co^5r)^(/>f ||ii|||.(j,3)||Vti|||.(^3)^^xrfs)i 

< (c|^sr)§(/>i5||K||i.(j,3)Cixd.)i(/>iV^||i.(^3)Cl^^^^ 

< (Co^Sr')lA 

The other integrals on the right side of (3.18) are bounded in a similar way, and (3.13) 
follows. 

For 1 < t < T, as in (3.18), we have 

sup /^3 \u\^ + \Vd\^dx + /ig3(|VMp + \V^d\^)dxds 

l<s<t 

< Mil! /m3 I Vci|>p(| V^xp + \V'd\') + m^Wdxds (3.19) 
+ /i /m3(I Vcipltip + I Vc^n Vc^tp + I Vdi|>|2)dxcls 
+Co + E]+A{l). 

The terms on the right can be bounded by 



/i /r3 \Vd\^\u\^\Vu\'^dxds 
< IiL:^i\^d\' + H' + \^Odxds 



10 , 14 



< E+ sup \\i%Vd)is)r\ j;\\iu,Vd){s)\\l^^^,^ds 

l<s<t L 3 ( ■ 



< £^ + ^sup[||(^x,Vci)(5)|||.(^3)||V(u,V6?)(s)||i.(^3)] 

X /; II {u, Vd) (s) ||i(«3) II V(«, Vd) {s) Wl^^s^ds 

< E + MCqAE. 
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By (1.1c) and (3.12), we have 

/i /k3 \'^d\^\dt\'^dxds 

< M J' J^^ilVdl^lV^dl^ + \Vd\''\u\^\Vd\^ + \Vd\^)dxds 

< M J^,{\Vd\*\V^d\^ + \V^d\^\u\^\Vdf)dxds 

< ME + MCoAE. 

The last term on the right side in (3.19) can be bounded by 

^ /i(/m3 \^d\^dx)^f^, \ufdx)lds 
< sup /jj3 \ Vd\^dx)3 Jjg3 \ Vu\'^dxds 

l<s<t 

2(9-3) _ 2 

The other integrals on the right side of (3.19) arc bounded in a similar way, and (3.14) 
follows. Thus the proof of lemma is completed. 

Next we derive a number of auxihary estimates needed to close the bounds in the 
previous two lemmas. We begin with a bound for the vorticity u. 

Lemma 3.5. Assume that the hypotheses and notations of Proposition 3.1 are in 
force. Then for < t < 1 A T, 



sup L3 (7^(1 VF|2 + \Vuj\'')dx + £ La (^(1 VF|2 + \Vuj\^)dxds 

0<s<t 

< M{D + eIc^^B^ + sup (T^\u\^dx + Jo (T\u\'^dxds), 

0<s<t 

and if T > 1 and 1 < t < T, then 

sup /k3(|VFP + \Vou\'')dx + J' /k3(|VF|2 ^ |Va;|2)(ix(is 



(3.20) 



< M{D + E + sup J^3\u\'^dx + J^3\u\'^dxds). 

l<s<t 



Proof. By (2.10) and the definition of D,E, we can easily get (3.20) and (3.21). The 
proof of lemma is completed. 

Next we derive an estimate for the functional Bg. 

Lemma 3.6. Assume that the hypotheses and notations of Proposition 3.1 are in 
force. Then for any < t < T, 

p — q 2 5+3 39—11 9—3 39 — 5 

Bg < M{C^-' + Bf^ + Cf^ Bf^ ) . (3.22) 
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Proof. We multiply (1.1b) by '^u and integrate over x (0,i) to obtain that 



/o /m3 - 2)|M|«-^((izt;w)xi • V{\u\^)dxds. 



For any > 0, 



/o /m3 |^(? - '2.)\u\'i-Hd;ivu)u ■ V{\u\'^)dxds\ 



< i%-2)/o7^3 



t r I J- Ml 1 



m| 2 2 |V(|'u| )|da;(is| 



< lX{q-2)[r,rj^s\ur'\{divu)\'dxds + r,-'j;j^s\ur'MW^^^^^ 
so if we choose 

^A(g-2)r7 = /3/x + A 
for a positive /3 to be determined, then the term in question will be bounded by 

3/3//, /q /jg3 \u\'^~^\Vu\^dxds + X /q /j^s \u\'^~'^\divu\^dxds 

+^-^t^£Ls\ur'mun'dxds. 

Substituting this into (3.23), we then get 

g-i ^3 p\u\'idx\l + p{l - 3/?) /J ^3 \u\i-^\Vu\^dxds 
+ [i/i(g - 2) + ti^^] ^3 l^il^-^l Vd^npdxd. 
< I /o /iR3(i^(p) - P(p))rf^^;(|«|'''2«)^^xrf5| 
+ 1 /o /m3 l«l^~^^Vrf • Addxds\. 

Recall that q G [Q,p), thus (1.9) holds with p replaced by g, and this is the condition 
that brackets on the left here is positive when /3 = |. It follows this term is positive 
for some /5 e (0, ^), which we now fix. It then follows that 



q-'J^,p\u\'^dx\l + J^f^, \ur'\Vu\'dxds + J^J^, \ur'M\u\')\'dxds 
< M[\ J^,{P{p) - P{p))div{\ur'u)dxds\ + I £ lu^'uVd ■ Addxds\]. 

(3.24) 

We multiply (3.9) by |Vo?|^"^V(i and integrate over x (0,^) to obtain that 

g-i /^3 \VdWi + £ ^3 iVdr'lV^-'dxds 
+ JiUq-2)\Vdr'\V{\Vd\')\'dxds (3.25) 
= /o /m3 \Vd\i-^WdWi\Wd\^d - u ■ Vd)dxds. 
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Adding (3.25) to (3.24) and applying the Cauchy's inequality in an elementary way we 
then obtain 

l^^l^ + \Vd\idx + £ \u\'^-^\Wu\^dxds + £ /j,3 \u\i-^\Vi\u\^)\^dxds 

+/o7m3 \^dr'\v'd\'dxds + /„7^3 \vdr^\vi\vdn'dxds 

< M[/j,3 Iwol" + \Vdo\^dx + I £ f^,{P{p) - P{p))dzv{\ur'u)dxds\ 

+ 1 /o /r3 \u\^~^uVd ■ Addxds\ + I \Vd\'i-^VdV{\Vd\^d - u ■ Vd)dxds\] 

(3.26) 

Since q e [6, min{p, 12}), then by Holder's inequality and Sobolev's inequality, we have 

/i < ( / \uq\^ + \VdQ\^dxY^^{ [ \uo\P+\Vdo\Pdx)^ < Mcl^'Ni^, 

and 

h < [Jq /«3 \u\'''-'dxds]m^ \Vu\'dxds]'^ 

< C|[/o(A3 l«P^c/x)l(/«3 WlKo-'Uxf^ds]'^ 

< C'|[/o(/m3 \ur'\Vu\'dx)iJ^, Hli^-^^dxy^ds]-^ 

1 i 12-q „ 39-16 

< C|5| sup [/j,3 |M(s)|2rfx]^ /„3 |^(s)|''da;]^ 

0<s<t 

9+3 3(j-ll 

< MC^B^. 

Using the fact 

Vrf • Arf = c/w(Vc/ (8) Wd) - ^Vdl Vd|^ 
and integrating by parts, we have 

^3 < [/o7^3(|«|'^-^|V«p+ |«|'^-^|V(|«p)P)^^x^i5]^[£43 |V(i|>|'^-2dxd.]^ 

< Mi?|[(£43 |Vd|^+2dxrf,s)l + (/*43 hvf+^dxd.)!] 

< ^4[(/o(/m3 |Vd|=^«dx)i(^3 |yd|='dx)ici.)^ 
+ (/o(/|3 l«l'j_f^a;)3(/^3 |Mpo?x)i(is)^] 

Similarly, we have 

h < M[/p7^3 \V'd\\Vd\'^-\\Vd\^ + \u\\Vd\')dxds] 

< M[/„*/^3(|Vrf|^^-2|V2rfPcixci.]^[/J/^3 \Vdr' + |V^in^.p^ia;ci5]i 

Substituting these results into (3.26) gives (3.22). Thus the proof of lemma is com- 
pleted. 

Next we derive a bound for the functional D and D. 
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Lemma 3.7. Assume that the hypotheses and notations of Proposition 3.1 are in 
force. Then for < t < 1 A T, 

D < M[C^B^A + A'^ + A^], (3.27) 
and if T > 1 and 1 < i < T, then 

D < M[cIA^ +A'^ + A% (3.28) 



Proof. We give the proof of (3.27), that of (3.28) being similar. By Lemma 2.1, we 
have 

< ^'l|Vci(-,t)||L[^l|VM-,t)lli.(M3)] 

< MA[a4||Vci(-,t)||i,(j,3) +a4||VM-,0lliW 



< MA[cpBr + lk^V2d(-,t)|||.(^^3)lkivM-,^)ll!.(M3)] 

< MA[Co*^S|^ + Ai||afvM-,^)lli(M3)]- 
Prom (3.9) and (3.12), we have 

||aiv3d(-,^)||i,(^3) < Mj^,a^{\Vdt\^ + \u\^\V^d\^ + \Vu\^\Vd\^ + \V^d\^\Vd\^)dx 
< M{A + D), 

so that 



/ 



a^\V''d\^\Vdfdx < MAlci"-' b/-' + A-^{A + D)-^] 



The other terms included in D are estimated in exactly the same way, and (3.27) fol- 
lows. The proof of lemma is completed. 

The following lemma contains the required bound for the pressure term in (2.11), 
which has been proved in Hoff [8, Lemma 3.3]. 

Lemma 3.8. Assume that the hypotheses and notations of Proposition 3.1 are in 
force. Then it holds 

/ / a'\p - pl^dxds < M[Co + [ [ a^\F\^dxds. (3.29) 
Jo Jut^ Jo Jr^ 

We can now obtain the required estimates for the functional E and E. 

Lemma 3.9. Assume that the hypotheses and notations of Proposition 3.1 are in 
force. Then there arc polynomials ipi and whose degrees and coefficients depend 
on the same M quantities as M in the statement of Proposition 3.1 such that: for 
0<t<lAT 

E < M[ip,{Co) + + B,)], (3.30) 
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and if T > 1 and 1 < i < T, then 

E < M[ipi{Co + A{1) + B,{1)) + <p2{A + Bg)]. (3.31) 

The polynomial cpi contains no constant term and the monomials in <^2 all have degrees 
strictly greater than 1. 

Proof. Since the term | Yl lo Ir^- ^'^xi "^^i '^xi, dxds\ has been bounded exactly 

1 )Jm ^3 

in Hoff [8]. So here we just bound the other terms for simplicity. 
First for < t < 1 A T, from (2.1), (3.9) and (3.12) we have 

/o Ir-^ cr^V'^d\^dxds 

< M(/* /jj3 a\V^d\dxds)-^{J^ ^3 a|V3rf|2rfxrfs)3 

< MA^J^ ^3 c^il^dtl^ + \Vd\^\V^d\^ + |«|2|V2rf|2 + \Vd\^\Vu\^)dxds)-^ 

< MA-*{A + /J(43 |Vd|6 + \uf + al{\S/''d\^ + \Vu\^)dxds)-* 

< MA{A + ct^Bf^ +E)-4, 

and 

J^J^,a^\V'd\^dxds 

< mlo^'iLs |V^dM:r)i(/^3 \V^'dx)lds 

< MiJ^{aJ^,\Vy\d^)HaJ^, \W^'dx)Ha' J^, \W^'dx)Us 

< MA{A + C^Bf^ + E)-4{A + D)-4. 

Prom Lemma 2.4, Lemma 3.2, Lemma 3.5 and the definition of F,uj, we have 

lo Ir3 cr^l^ul^dxds 

< M[/* ^3 a'{\p -~P\' + \F\'+ \u\^)dxds] 

< M[Co + (sup J^,a{\F\' + \u;\')dxJ^,a%\VF\' + \Vu\')dx)h 

0<s<t 

x{J^J^s<T{\VF\' + \Vu;\')dx) 

< M[Co + (Co + A)-2{A + D)-2{A + Elcf^' B^')], 

lo lR3^h^^\^dxds 

< M[J^J^,aH\p-p\^ + \Ff+\uf)dxds] 

< M[C'o + />t(^3(|FP + MW^(/M3(|VFP+|Va;nda;)id. 

< M[Co + {Co + A)l{A + Eicf^'B^')l]. 

Thus combining the above results and Lemma (3.7), we yields (3.30). 
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Now for 1 < f < T, if we take g = 4 in (3.25) and integrate by parts to obtain that 

^3 \Vd\^dx + /; ^3 \Vd\^\V^d\Hxds + jl ^3 \V{\Vdn^dxds 

< M[/j,3 I W(-, + I jl ^3 l\d\Vd\^ddxds\ + jl ^3 \Vd\^\V^d\ \u\dxds\ 

< M[C(r^Sg(l)^ + /^3 1 Ad + \Vd\^d\^\Vd\^dxds 
+(/o 4' \^d\f\V^d\5dxds)i{Jl^J^, \ufdxds)^. 

< M[Co^5,(l)i^ + sup II Vrf||i^.jj3) /* L I Ad + I Vdpdpdxds 

l<s<t 

+(/o 1 vdpi v^dpdxd.)! J5,3 1 vtxpdxds)^ 

Multiplying (3.9) by VAd and integrating over R^, we have 

\V^^dx < ^3 \Vdt\\V^dx + 2 ^3 \Vd\\V'd\\V^dx 
+ /jj3 Vm ■ VdVAd + w ■ WdVAddx 
= /e3 |V(it||V3d|(ix + 2/^3 |Vd||V2(i||V3(i|dx 

+ Vu ■ VdVAd - Vu ■ WdAd+ l{divu)\Ad\'^dx. 

By Cauchy's inequality, we have 

iV^-'dx < /j,3 iVdtP + iVdPlV^c/p + |Vdp|V«P + iViillV^d^dx. 

Thus we have 

/i /e3 iV^dl^dxds 

< M(/;(/^3 |V2d|dx)i(/^3 iV^dpda:)!^. 

< MA/^*(/j,3 |Vdt|2 + iVdplV^dp + |Vd|2|VMp + \Vu\\V^^dx) 
x(/r3 iVdip + iVdHV^dp + iVdplVwp + |^i|2|V2(i|2(ia;)5ds 

< MA{A + 5)^(/i (/^s iVdtP + \Vd\''\V^''dxds 
+ f.iLs \Vd\'\Vu\'+^u\\V'd\'dxds) 

< MA{A + D)-2{A + cl^Bq{l)^^ 
+Co(C*^S|^ +Ai{A + D)i) + cIeI + cIe-2). 

Bounds for the term /,* /jjadV-up + |Vii|^)da; are obtained in a similar way, which in 
fact is much more simple. Then applying Lemma 3.7, we can bound E which gives 
(3.31). The proof of lemma is completed. 

Combining the results of Lemmas 3.2-3.9, we have the following bound for A + Bg. 

Lemma 3.10. Assume that the hypotheses and notations of Proposition 3.1 are in 
force. Then there are polynomials ifi and y?2 as described in Lemma 3.9 such that for 
< t < 1 A T, 

A + Bg< M[<p^{Co) + MA + B^)], (3.32) 
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and if T > 1 and 1 < i < T, then 

A + Bg< M[ip,{Co + A{1) + B,{1)) + if2{A + B,)\. (3.33) 



Proof of Proposition 3.1. Proposition now follows immediately from the bounds 
(3.32) and (3.33) and the fact that the functions A, A, Bq are continuous in time. 

4 Pointwise bounds for the density 

In this section we derive pointwise bounds for the density p, bounds which are 
independent both of time and of initial smoothness. This will then close the estimates 
of Proposition 2.1 to give an uncontingent estimate for the functional A defined in 
(3.1). 

We begin with two auxiliary lemmas. The first lemma is a maximum-principle ar- 
guments applied integral curves of the velocity field, which has been proved in Hoff [8]. 

Lemma 4.1 Let {p,u,d) be as in Proposition 3.1 and suppose that < Ci < p < C2 
on X [0, T]. Fix to > and define the particle trajectories x : [0, oo) x ^ by 

xit.y) = u{x{t,y),t), 
x{to,y) = y. 

Then there is a constant C depending only on Ci and C2 such that if 5* G L^(]R'^) is 
nonnegativc and t G [0, T], then each of the integrals J^^ g{x{t, y))dy and Jj^g g{x)dx is 
bounded by C times the other. 

Next we derive a result relating the Holder-continuity of u{-,t) to various norms 
appearing in the definition (2.1) of the functional A. 

Lemma 4.2. Let {p,u,d) be as in Proposition 3.1. Then for a G (0, |] and t G (0,r], 
we have 

{U{;t)r < M[||Vtx(-,t)||^3)||Vc.(.,t)||^3) + (Co + \\Vui; 

xiM-Mhm + \\i^d-Ad){;t)\\hmy-^ + 



Proof. Let a G (0, i] and define r G (3,6] by r = ^z^. Then by (2.3) and (2.11), we 
have 

{ui;t)r < M[||F(-,t)|U.(M3) + ||a;(-,0|U.(M3) + ||(p - p)(., t)||^.(j,3)]. (4.2) 
By (2.1), we obtain 

M;t)\\Lrm < ^(IM->i)llS(iR3)l|Vc^(-,t)||5"f 



< ^l|V«(-,t)||^.V3)||Va;(-,t)||^4^3)) 
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and 

||i^(-,i)||L^(M3) < M(||F(-,t)|||p,3)||VF(-,t)||^,(^3); 

< M(||(p-p)(-,f)||2 +||V<,f)||i. 



x(||ii(-,t)|li.(R3) + ||(Vd-Aci)(-,t)||2,(^3))^ 
< M(||Co+||V<,t)||i.(«3))^ 
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x(IK,t)|li.(M3) + ll(Vrf-Aci)(-,t)||i.(^3))- 

Putting the above results into (4.2) yields (4.1). Thus the proof the lemmas is com- 
pleted. 

Now we derive the upper and lower pointwise bounds for the density. 

Proposition 4.3. Assume that the system parameters in (1.1) satisfy the condi- 
tions (1.5)-(1.8) and let positive numbers and 6 < 5 be given. Assume [p,u,d) 
is a solution of (1.1) on x [0,T] in the sense of Proposition 2.3 with initial data 
(pc^o) e H%R^) and do E H^{R^; E>^) satisfying (1.10)-(1.13). Then there are positive 
constants e, M, and 9 depending on the parameters and assumptions in (1.5)-(1.8), on 
N, and a positive lower bound for b, such that, if Cq < £ and p{x, t) > on R3x[0,T], 
then in fact 

P<P<P on X [0,T], (4.3) 

and 

A{T) < MCl (4.4) 



Proof. First we choose positive numbers k and k' satisfying 

p<n<p + b<p — b<K'<p. 

Recall that po takes values in [p + 6, p — 6], so that p G [p, p] on M? x [0,r] for some 
positive r by the time regularity (2.7). It then follows from Proposition 3.1 that 
A{t) < MCq, where M is now fixed. We shall that if Cq is further restricted, then in 
fact that K< p< k' on all of x [0, T], and therefore that A{T) < MCl as well. We 
shall prove the required upper bound, the proof of the lower bound being similar. 
For y e and define the corresponding particle path x{t) by 

f x{t,y) = u(x(t,y),t), 
1 x{to,y) = y. 

Suppose that there is a time h < t such that p{x{ti),ti) — k'. We may take ti minimal 
and then choose to < ti maximal such that p{x{to),to) — p — b. Thus p{x{t),t) G 
[p — b, K.'] for t E [to, ^i]- We divide into two steps: 
Step 1. to<ti<TAl 

We have from the definition (1.3) of F and the mass equation that 

+ X)^\\ogp(x{t),t) - \ogp] + P{p(x(t),t)) - P(p) = -F{x{t),t). 
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Integrating from to to ti and abbreviating p{x{t),t) by p{t), etc., we then obtain 

{p + \) log p{s)\ll+ f\p{s)-P{p)]ds = - r F{s)ds. (4.5) 

J to J to 



f 

J to 



We shall show that 

F(s)ds < MCl (4.6) 

'to 

for a constant M which depends on the same quantities as the M from Proposition 3.1 
(which has been fixed). If so, then from (3.5), we have 

{p. + A) [log k' - log(p -h)]<- r [P{s) - P{p)]ds + MCl < MCl (4.7) 

Jto 

where the last inequality holds because p{t) takes values in [p — b, k'] C [max{p, p"}, p], 
and P is increasing on [p, p\. But (4.7) cannot holds if Cq is small depending on M, k', 
and p — b. Stipulating the smallness condition, we therefore conclude that there is no 
time ti such that p{ti) — p{x{ti),ti) = k' . Since y was arbitrary, it follows that 
p < on X [0, r], as claimed. The proof that p > k is similar. 

To prove (4.6) we let F be the fundamental solution of the Laplace operator in 
and apply (1.4) to write 

Ko = J^,{V^T{x{s) - y))puiy, s)dyds 

+ mA^.nx{s)-y)){Vd-Ad){y,s)dyds. ^ ' ' 

By Lemma 4.2, the first integral on the right side of (4.8) is bounded exactly as in 
Lemma 4.2 of Hoff [8]: 

Ito /R3(Va;r(x(s) - y))pu{y, s)dyds 

< ||Vr * (pti)(-,tl)||L°°(K3) + II Vr * (pM)(-,t2)||L-'(R3) 

+ /o /k3 ^x,xk{x{s) - y)[u^{x{s), s) - u''{y, s)]{pu^){y, s)dyds 

< MC'o + MC'oJ^{u{-,s))^ds 

< MCo^ + MCo^(^S-^ds)i/^^(C'o+ ||V«(,t)||i,(^3)rf^^ 
x(/o^nH-,t)|li.(M3) + ll(Vci-Ad)(-,t)||i,(^3))rf.)^ 

< MCo^ 

if a < |. Note that (3.22) holds for q = 6, thus if 2 < r < by (2.4), the second 
integral on the right side of (4.8) can be bounded in as 

l!o /R3(V.r(x(s) - y))iVd ■ Ad)iy, s)dyds 

< MJ^ ||(VdAd)(s)|U2(K3) + ||(VrfArf)(s)|U.(K3)rfs. 

< MJ^ ||(|VdnArfp)(s)||i,(^3)|||Adp(s)||i,(^3)C^s 

+ /O II ( Ad) (S) |U3(M3) II ( Vd) (S) 11^^ ^^3^d. 
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Thus the proof of (4.6) is completed. 
Step 2. 1 < to < ^1- 

Again by the mass equation and the definition (1.3) of F, 

j^iPit) - p) + (/. + xr'pitmt) -p) = {p+ xr'p{t)F{t). 

Multiplying by {p{t) — p) we get 

\{p{t) - pf + (/. + \r'f{t)p{t){p{t) - pf = -{p + X)-'p{t){p{t) - p)F{t), (4.9) 
where 

m = (Pit) - p)ipit) - p)-\ 

Since f{t) > on [tc^i], thus integrating (4.9) over [toj^i], we arrive at 

\^'-p\^-\p-b-p\'<M \\F{; s)\\'s.ds. (4.10) 



to 



So that if we show that 



r \\F{;s)\\l^ds<MCl (4.11) 

J to 



Then as in Step 1, (4.10) cannot hold if Co is sufficiently small. Since ^/ e was 
arbitrary, it follows that p < /t' on x [0,r], as claimed. 
To prove (4.11) we apply (1.4) and (2.4) to get 

Ko \\F{;s)\\l^^^.^ds < \\u{-,s)\\l,^^,^ds + !ll IKVrf- Ad)(-,.)||i,(^3)C/. 

ii^(->^)iii4(M3)rf^+/;; \\iyd-Ad){.,s)\\i,^^,^ds 
< Mc^, + \\u{; .)iii.(^3)rf. + iKvrf • Ac/)(-, .)ii 



The terms integral on the right side above can be bounded as 

ds < ||'!i(-,s)|l!2(K3)||Vlt(-,s)|||2(M3)C?S 

< ui: \\u{:s)\\i,^^,^ds)Hsii iiv^(-,.)iii.(^3)rf.)t 

< MC^Q 

and 

/^MI(Vd-Ac/)(-,5)||i,(^3)d. 

< ||(vdAd)(-,.)|ii,(^3)||v(vd. Ad)(.,.)||f,(^3)rfs 

< {Jll\\{vd-/:^d{-,s)\\l,^^,^ds)l{jl^j^^ 

where the last inequality follows from Proposition 3.1. Thus (4.11) is proved. The 
proof of Proposition is completed. 
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5 Proof of Theorem 1.2 



In this section, we prove Theorem 1.2 by constructing weak solutions as limits of 
smooth solutions. So, we first prove the global-in-time existence of smooth solutions 
with smooth initial data which is strictly away vacuum and is only of small energy. 

Proposition 5.1. Assume that {po,Uo,do) satisfy (2.5). Then for any < T < oo, 
there exists a unique smooth solution {p,u,d) of (1.1)-(1.13) on x [0,r] satisfying 
(2.6)-(2.9) with Tq being replaced by T, provided the initial energy Co satisfies the 
smallness condition (1.16) with £ > being the same one as in Proposition 3.1 and 
Proposition 4.3. 

Proof. The standard local existence result (Proposition 2.3) shows that the Cauchy 
problem (1.1)-(1.2) admits a unique local smooth solution {p,u,d) on IR^ x [0,To]. In 
view of Lemma 3.2 and Proposition 4.3, we have 

A{To)+ sup J^,i\p~p\^ + \u\' + \Vd\^)dx 

0<t<TQ (5.1) 

+ /o^"/r3(|VmP + |A(i+ |Vd|2(i|2)M < MCo, 

and 

P< P<P on X [0,To]. (5.2) 

Then the standard arguments based on the local existence results together with the a 
priori bounds (5.1)-(5.2), we deduce that (p, u, d) is in fact the unique smooth solution 
of (1.1)-(1.13) on X [0, T] for any < T < oo. 

With the help of Proposition 5.1, we are in a position to prove Theorem 1.2. 

Proof of Theorem 1.2. For any map do e Hl{'^^]'S^), there exists d'^' € H^{R^-'^'^) 
such that 

lim II - c?o||iji(]R3) = 0. 

m— >-oo 

Let 

m m 

where Jj_ = Jj_(x) is the standard moUifier. Then (p[J^,w^) G i/^(M^) and d"^ — n E 

m m 

//^(M^;§2) and the initial norm for {p^,ul^,Vdl^) (i.e., the right side of (1.13) with 
(po, 1*0) Vo?o) replaced by (p^, u^, Vd^)) is bounded by Cq. The above Proposition can 
be applied to obtain a global smooth solution {p"^ , u"^ , d^) of (1.1)-(1.13) satisfying 
(3.2), (4.3) and (4.4) for alH > uniformly in m. 

In view of (2.3) and (2.11), we sec from Sobolev embedding theorem that 

(w"^(-,t))^ < C||Vii™|Ua(M3) 

< C(||F-|U.(K3) + ||a;'"||L6(R3) + IIP- - P|U6(M3)) 3^ 

< C(r)(l+||M™|U2(K3) + ||Vd-- Ad™|U2(K3)) ^'^ 

< C{t), t>T>0, 
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where F"*, a;"* and are the functions F, a; and P with (p, u, d) being replaced by 
In addition to (5.3), we also have 



u"^{x, t) - f u"^{y, t)dy\ < C{t)R\ 



and hence, for < r < ti < ^2 

< \Bk:;\\f:!B,Jn"'iy^t)\dydt + c{r)R\ 

< CR-i\t, - t,\Hj;; J^^^^^ \u"^{y, t)\^dydt)h + C{t)R^ (5.4) 

< cR-h2-h\Hp;j^Jj\ir{y,t)\' + \u^\'\wu\')^^^^ 

< C{T)[R-l\t2-ti\-2 + R^. 
Taking i? = — in (5.4), we get 

\u"'{x,t2)-u"'{x,ti)\<C{T)\t2-ti\K 0<T<ti<t2<OO. (5.5) 

The same estimates in (5.3) and (5.5) also hold for d and Vrf. Thus, wc have proved 
that {w™}, {c?™} and {Vrf™} are uniform Holder continuity away from t = 0. As a 
result, it follows from Ascoli-Arzela theorem that 

u"^ — > u, d^ — V d uniformly on compact sets in M.^ x (0, oo). (5.6) 

Moreover, by argument in [24] (see also [5]), we know that 

— > p strongly in LP{R^ x (0,oo)), Vp e [2,oo). (5.7) 

Therefore, passing to the limit as m — > oo by (5.6) and (5.7) we obtain the limited 
(p,u,d) which is indeed a weak solution of (1. !)-(!. 13) in the sense of Definition 1.1 
and satisfies (1.17)-(1.24). 

Next we derive the large-time behavior of {p,u,d) in (1.25). This can be done as 
the ones in [4] , however, for completeness we sketch the proof here. We first deduce 
from the mass equation that 

(P(p) -P)t + u- V(P(p) - P) + -iP{p)divu = 0. 

Multiplying the above equation by 4:{P{p) —P)^ and integrating it over M^, we get that 

\\P{p) - P\\Um = I mp) - P\')divu-3^P{p){P{p) - Pfdtvudx, 



7,11- vr/ - ui.-(iKi) 

which, together (3.29) shows that 



/oo ^ poo POO 

\-\\P{p)-P\\l,^^s)\dt<C{l + \F\l,^^s^ds) |Vw|i2(ij3)fl!s < C, 
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As a result, we have 

\\P{p) - ^||l4(m3) — >0 as t^oo. 

This, together with (3.2) and the uniform lower and upper bound of density, shows 
that 

lim ||p-p||ii(R3) = (5.8) 



holds for any / G (2, oo). 

Following the argument in [4], wc take a sequence 

u"'{t,x) := u{t + m,x), 

for all integer m, and {x,t) e x [1, 2]. Then from (1.24), we have 



-1 

lim / \Wu"'\\r2ms\ = 0. 







Prom (1.24) again, we have 

||'f^"*||ifi(R3) < C uniformly for t,m. 

Thus we arrive at 

lim ||M"^||L2(]g3) = uniformly for t. 

That means 

lim ||«(i)||i2(M3) =0. (5.9) 
For t > 1, from (2.10) and (2.11), we obtain that 

||VM(t)||L6(K3) < C(||F(t)|U6(R3)) + ||a;(t)|U6(K3) + ||(P(p) - P))(0||l6(R3) 

< C(l + ||VF(t)|U2(M3) + ||Va;(t)|U2(K3)) . . 

< C(l+||«(t)|U2(M3) + ||(VdAd)(t)|U2(K3)) ^ 

< c. 

Combining (1.24) (5.9) and (5.10), we have 

lim ||M||vyi.'-(R3) = 0, (5-11) 



t— >-oo 

holds for r G (2, 6). 
Similarly, we have 

lim ||Vrf||H^i..(M3) = 0, (5.12) 

t^oo 

holds for r G (2,6). Putting (5.8) ,(5.11) and (5.12) together gives (1.25). Thus the 
proof of Theorem 1.2 is completed. 
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